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We suggest to use insulating garnets doped by rare earth ions for measurements of nuclear anapole 
moments. A parity violating shift of the NMR frequency arises due to the combined effect of the 
■ lattice crystal field and the anapole moment of the rare-earth nucleus. 

We show that there are two different observable effects related to frequency: 1) A shift of the 
NMR frequency in an external electric field applied to the solid. The value of the shift is about 
A^i ~ 1CT 5 Hz with E = 10 kV/cm; 2) A splitting of the NMR line into two lines. The second 
effect is independent of the external electric field. The value of the splitting is about A^2 ~ 0.5 Hz 
and it depends on the orientation of the crystal with respect to magnetic field. Both estimates are 
presented for a magnetic field of about 10 tesla. 

We also discuss a radiofrequency electric field and a static macroscopic magnetization caused by 
the nuclear anapole moment. 
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I. INTRODUCTION 



The anapole moment is a characteristic of a system which is related to the toroidal magnetic field confined within the 
system. It was pointed out some time ago by Zeldovich Q that the anapole moment is related to parity violation inside 
' the system. Interest in the nuclear anapole moment is mostly due to the fact that it gives dominating contribution 
to effects of atomic parity nonconservation (PNC) which depend on nuclear spin 0. There are two mechanisms 
that contribute to these effects. The first is due to exchange of a Z-boson between electron and nucleus. The 
j^n V second mechanism is due to the usual magnetic interaction of an electron with the nuclear anapole moment. The 
contribution of the first mechanism is proportional to 1 — 4s 2 . Since sine squared of the Weinberg angle is s 2 « 0.23 
0, the first mechanism is strongly suppressed and the second mechanism dominates. The anapole moment of 133 Cs 
has been measured in an optical PNC experiment with atomic Cs 4] . This is the only observation of a nuclear anapole 
moment. There have been several different suggestions for measurements of nuclear anapole moments. Measurements 
in optical transitions in atoms or in diatomic molecules remains an option, for a review see [{|. Another possibility 
■^■j- ■ is related to radiofrequency (RF) transitions in atoms or diatomic molecul es pi , 0, H, 0- Possibilities to detect 
nuclear anapole moments using collective quantum effects in superconductors jlOj . as well as PNC electric current in 
ferromagnets |Tlj , have been also discussed in the literature. A very interesting idea to use Cs atoms trapped in solid 
4 He has been recently suggested in Ref. [l2|. 

Our interest in the problem of the nuclear anapole moment in solids was stimulated by the recent suggestion for 
searches of electron electric dipole moment in rare earth garnets |l3j . Garnets are very good insulators which can 



^ , be doped by rare earth ions. They are widely used for lasers and their optical and crystal properties are very well 



understood. To be specific we consider two cases: the first is yttrium aluminium garnet (YAG) doped by Tm [1^ 



Thulium 3+ ions substitute for yttrium 3+ ions. The second case is yttrium gallium garnet doped by Pr 15]. Once 
more, praseodymium 3+ ions substitute for yttrium 3+ ions. The dopant ions have an uncompensated electron spin 



J and a nuclear spin I. For Tm 3 + J = 6 and 7=1/2 ( 169 Tm, 100% abundance). For Pr 3+ J — A and I = 5/2 
( 141 Pr, 100% abundance). 

The simplest P-odd and T-even correlation (P is space inversion and T is time reflection) which arises due to the 
nuclear anapole moment is 
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where E is the external electric field. It is convenient to use the magnitude of the effect expected in the electron 
electric dipole moment (EDM) experiment [HI as a reference point. For this reference point we use a value of the 
electron EDM equal to the present experimental limit 0], d e = 1.6 x 10~ 27 e cm. According to our calculations, the 
value of the effective interaction is such that at the maximum possible value of the cross product [I x J] it induces 
an electric field four orders of magnitude higher than the electric field expected in the EDM experiment 0, . 
For example, in P^CasO^ the field is E ~ 1.5 x 10~ 6 V/cm. The problem is how to provide the maximum cross 
product [J x J]. Value of (J) is proportional to the external magnetic field B. A magnetic field of about 5-10 T is 
sufficient to induce the maximum magnetization. Nuclear spins can be polarized in the perpendicular direction by 
an RF pulse, but then they will precess around the magnetic field with a frequency of about 1 GHz. It is not clear 
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if the anapole-induced voltage of this frequency can be detected. An alternative possibility is to detect the static 
variation of the perpendicular magnetization induced by the external electric field, SI oc [B x E] . The magnetization 
effect for Pr3Ga50i2 is several times larger than that expected for the EDM experiment 0]. This probably makes 
the magnetization effect rather promising. In the present work we concentrate on the other possibility which is based 
on the crystal field of the lattice. Because of the crystal field, the electron polarization of the rare earth ion has a 
component orthogonal to the magnetic field (J) oc B + (B ■ n)n, where n is some vector related to the lattice. The 
equilibrium orientation of the nuclear spin is determined by the direct action of the magnetic field together with the 
hyperfine interaction proportional to (J). Because of the (B ■ n)n term in (J), the nuclear and the electron spins are 
not collinear, and the cross product [I x J] is nonzero [I x J] oc (B ■ n)[B x n]. We found that NMR frequency shift 
due to the correlation is about 

Avi ~ 1(T 5 Hz (2) 

at E — 10 kV/cm and B = 10 T. In essence, we are talking about the correlation (B ■ n)[B x n] ■ E considered 
previously in the work of Bouchiat and Bouchiat 0] f° r Cs trapped in solid 4 He . 

Another effect considered in the present work is the splitting of the NMR line into two lines due to the nuclear 
anapolc moment. This effect is related to the lattice structure and is independent of the external electric field. 

The garnet lattice has a center of inversion. However, the environment of each rare earth ion is asymmetric 
with respect to inversion. One can imagine that there is a microscopic helix around each ion. Since the lattice is 
centrosymmetric, each unit cell has equal numbers of rare earth ions surrounded by right and left helices (there are 24 
rare earth sites within the cell). The microscopic helix is characterized by a third rank tensor Tkim (lattice octupole). 
Together with the nuclear anapole interaction this gives a correlation similar to , but the effective "electric field" 
is generated now by the helix oc T^i m JiJ m . So the effective interaction is 

H^g OC eijkhJjTklmJlJm- (3) 

The effective interaction @ produces a shift of the NMR line. The value of the shift is about 0.5 Hz at B = 10 T, 
and the sign of the shift is opposite for sites of different "helicity" , so in the end it gives a splitting of the NMR line 

Az/ 2 ~ 0.5 Hz. (4) 

The value of the splitting depends on the orientation of the crystal with respect to the magnetic field. This is the 
"handle" which allows one to vary the effect. Generically this effect is similar to the PNC energy shift in helical 
molecules fl8| . 

One can easily relate the values of the frequency shift in the external field (J2J and of the line splitting Q. The 
splitting is due to internal atomic electric field which is about 10 9 V/cm. Therefore, naturally, it is about 5 orders of 
magnitude larger than the shift (J2J) in field 10 kV/cm. 

For the present calculations we use the jelly model suggested in Ref. [IT]]. Values of the nuclear anapole moments 
of 169 Tm and 141 Pr which we use in the present paper have been calculated separately 0|. The structure of the 
present paper is as follows. In Section II the crystal structure of the compounds under consideration is discussed. The 
effective potential method used in our electronic structure calculations is explained in Section III. The most important 
parts of the work which contain the calculations of the effective Hamiltonians and are presented in Sections 
IV and V. The crystal field and the angle between the nuclear and the electron spin is considered in Section VI. In 
section VII we calculate values of observable effects and section VIII presents our conclusions. Some technical details 
concerning the numerical solution of the equations for electron wave functions are presented in Appendix. 



II. CRYSTAL STRUCTURE OF Y(Pr)GG AND Y(Tm)AG 

The compounds under consideration are ionic crystals consisting of Y 3+ , 2 ~, Ga 3+ ions for YGG and Al 3+ instead 
of Ga for YAG, plus Pr 3+ or Tm 3+ rare-earth doping ions. The chemical formula of YGG is Y^G^O^ and the formula 
of YAG is Y3AI5O12. Yttrium gallium garnet and yttrium aluminium garnet belong to the laid space group and 
contain 8 formula units per unit cell. Detailed structural data for these compounds are presented in Table [I] }20L |2l| . 

RE 3+ doping ions replace Y 3+ ions and hence enter the garnet structure in the dodecahedral 24c sites with the local 
D2 symmetry. In this case each RE 3+ ion is surrounded by eight oxygen O 2- ions in the dodecahedron configuration 
resembling a distorted cube (see Fig. ^) . There are 24 such sites per unit cell: half of them have absolutely identical 
environment with the other half; the remaining 12 can be divided into 6 pairs where the sites differ only by inversion, 
and these 6 pairs differ with each other by finite rotations. In the present paper we perform calculations for the 
case of one particular site orientation; the coordinates of the oxygen atoms around the central impurity ion for that 
instance are presented in Table |Hj After that, the results for all other sites in the unit cell can be found by applying 
the inversion of coordinates or the necessary rotations, listed in Table ITTll 
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YGG YAG 

unit cell parameters (A) 
a,b,c 12.280 12.280 12.280 a,b,c 12.008 12.008 12.008 

a,/?, 7 90° 90° 90° a, 0,7 90° 90° 90° 
space group 

laid (230 setting 1) laid (230 setting 1) 

atomic positions 
Y 0.1250 0.0000 0.2500 Y 0.1250 0.0000 0.2500 
Ga 0.0000 0.0000 0.0000 Al 0.0000 0.0000 0.0000 
Ga 0.3750 0.0000 0.2500 Al 0.3750 0.0000 0.2500 
O 0.0272 0.0558 0.6501 O 0.9701 0.0506 0.1488 



TABLE I: Structural data for YGG [U and YAG HH. 



FIG. 1: Dodecahedron configuration of O 2 ions around the RE 3+ impurity ion in the garnet structure. Two different viewing 
angles are shown. 



III. CALCULATION OF THE ELECTRONIC STRUCTURE OF RE O s CLUSTER 

We describe an isolated impurity ion with the effective potential in the following parametric form: 

r (l + w ) 2 {e^ +1) r 
Pr : fi= 1.0, d = 1.3, T] = 2.25; 
Tm : fi= 1.0, d = 1.0, r) = 2.56. 

Here Z is the nuclear charge of the impurity ion, Zi is the charge of the electron core of ion, and fi, d and rj are 
parameters that describe the core. We use atomic units, expressing energy in units of Eq = 27.2 eV and distance 
in units of the Bohr radius as — 0.53 x 10~ 8 cm. Solution of the Dirac equation with the potential 101 gives wave 
functions and energies of the single-electron states. The potential (J5J provides a good fit to the experimental energy 
levels of isolated impurity ions |23; the comparison is presented in Table llVl 
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TABLE II: Coordinates of oxygen ions in YGG and YAG (A) with respect to the rare earth ion. The axes x, y, and z are 
directed along the three orthogonal cube edges a, b, and c, Table □ 
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TABLE III: Euler angles of rotation between inequivalent RE 3+ impurity sites. 
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TABLE IV: Calculated and experimental |2^l energy levels of an isolated ion with respect to the ionization limit. Energy levels 
are averaged over the fine structure. Units 10 3 cm -1 . 

In order to model the electronic structure of the RE-Og cluster (Fig. QJ, following 17] we use the jelly model and 
smear the 8 oxygen ions over a spherical shell around the rare earth ion. Hence, the effective potential due to the 
oxygen ions at the RE 3+ site is 

Voir) = -A e-(^ 2 , (6) 

where r a = 4.5 as is the mean RE-0 distance, A a and D are parameters of the effective potential. To describe the 
electrons which contribute to the effect we use the combined spherically symmetric potential 

V{r)=V RE (r)+Vo(r), (7) 

where Vre is the single impurity ion potential Solution of the Dirac equation with potential JjJ gives the single- 
particle orbitals. In this picture we describe the electronic configuration of the cluster as [RE 3+ ]6s 2 6p 6 , where the 
electronic configuration of Pr 3+ is Is 2 ... 5s 2 5p 6 4/ 2 and Tm 3+ is Is 2 ... 5s 2 5p 6 4/ 12 . The eight states 6s 2 6p 6 represent 
2p CT -electrons of oxygens combined to S- and P- waves with respect to the central impurity ion (see Ref. |l7jh Parame- 
ter A in the "oxygen" potential Vo I© is determined by matching the wavefunction of oxygen 2p CT -orbital (calculated 
in Ref. j^) with the 6s- and 6p-orbitals from the combined potential Q at the radius R w 2.5as. The matching 
conditions are 

\MR)\ = \i>2 P Aro-R, cos e = i)\, 

\ipe p {R,cos6 = 1/V3)| = IV^r-o - -R, cos0 = (8) 

This is a formulation of the idea of dual description at r w R, see Refs. [2^112^ . 

The parameter D in © represents the size of the oxygen core and is about D S, 1 (atomic units). The jelly model 
is rather crude and the value of D cannot be determined precisely, see Ref. |l7| . In the present work we vary this 
parameter in the range of 0.1-1.5 . For each particular value of D we find A a to satisfy (JSJ), for example A Q = 0.9 at 
D = 1. The most realistic value for D is probably around 0.5-1.0. To be specific, in the final answers we present 
results at D = 1.0. Instead of the jelly model it would certainly be better to use a relativistic quantum chemistry 
Hartree-Fock method (or the Kohn-Sham form of the relativistic density functional method which allows one to 
generate electron orbitals) to describe the RE-Os cluster. However, this would be a much more involved calculation 
at the edge of present computational capabilities and therefore, at this stage, we continue with the jelly model. 

IV. CALCULATION OF THE EFFECTIVE HAMILTONIAN © 

The calculations in the present section are similar to those performed in |l7j for the electric dipole moment of the 
electron. There are three perturbation operators that contribute to the correlation ft}. First, there is a magnetic 
interaction of the electron with the nuclear anapole moment, see, e.g., |5j. Expressed in atomic units the interaction 
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reads 



V a 
K a 

141 Pr : n a 
169 Tm : k„ 



K a {Ia)6(r) 
, Jm 2 a 



(9) 



V2 



1.57-10- 14 K a S a , 



ID mro 

0.35, S'a = -0.34, 
0.39, S a = -0.25. 



Here m is the electron mass, G is the Fermi constant, and a is the fine structure constant; a are the Dirac matrices, \x 
is the magnetic moment of the unpaired nuclcon (proton in these cases) expressed in nuclear magnetons, fo — 1.2 fm, 
A is the mass number of the nucleus, and g « 4 for outer proton and g ~ 1 for outer neutron. Values of the nuclear 
structure constant S a have been calculated in |19j. 





FIG. 2: Schematic picture, illustrating the shift of Vo(r) due to the lattice deformation. 

The second perturbation operator is related to the shift Ar of the rare earth ion with respect to the surrounding 
oxygen ions. The shift is proportional to the external electric field, but for now we consider Ar as an independent 
variable. In the jelly model — Ar is the shift of the spherically symmetric oxygen potential Vfo(r) © with respect to 
the origin, see Fig. [21 Therefore, 



Vo(r)^V (r) = V (r + Ar) = V (r) + 

Thus, the perturbation operator related to the lattice deformation reads 

(Ar • r) dV Q 



(Ar • r) dV 

r dr 



Vi(r) 



dr 



= (Ax sin 9 cos 4> + Ay sin 9 sin 4> + Az cos 9) (—2 



(r- 



D 2 



(10) 



(11) 



Here r = r(sin 9 cos <fc, sin 9 sin <fi,cos9). 

The third perturbation is the residual electron-electron Coulomb interaction, which is not included in the effective 
potential, 



V c (ri,rj) 



1 



■in 



^21 + 1 r l + 



T F4(r,)Yi m (r i ). 



(12) 



Here r^ and Tj are radius-vectors of the two interacting electrons. 

The formula for the energy correction in the third order of perturbation theory reads, see, e.g., Ref. |26j |: 



£i 3) = 



n u 2 



(13) 



7) I 



where V = V a + V\ + Vc^ In Eq. we need to consider only the terms that contain all the operators V , Vi, and 

v c . 

The shift operator V\ is nearly saturated by 6s- and 6p-states because core electrons do not "see" the deformation 
of the lattice, hence, for this operator we consider only s-p mixing. Matrix elements of the anapole operator V a 
practically vanish for the electron states with high angular momentum, since this operator is proportional to the 
Dirac delta function. Therefore, it is sufficient to take into account only (nsi^l Vi lfcpi/2) matrix elements. All in all, 
there are 11 diagrams (Fig. that correspond to Eq. I|13|) . All diagrams are exchange ones and contribute with the 
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FIG. 3: Third order perturbation theory diagrams corresponding to Eq. 11311 . The cross denotes the anapole interaction V a ©, 
the dashed line denotes the lattice deformation perturbation Vi dill) , and the wavy line denotes the Coulomb interaction Vc 
1121 1. The multipolarity of the Coulomb interaction is shown near the wavy line. Each diagram contributes with the coefficient 
shown before the diagram (number of diagrams of this kind). Summation over all intermediate states |fc) and |m) and over all 
filled states \n) is assumed. 



sign shown before each of the diagrams. Summation over all intermediate states and \m) and over all filled states 
\n) is assumed. 

Since V a and V% are single-particle operators, we evaluate each diagram by solving equations for the corresponding 
wavefunction corrections. For example, the first diagram contains in the top right leg the correction 

^ x )=Y. {mPl ' MnS) \rn Pl/2 ). (14) 
m £ ™ s £m Pi/2 

To evaluate the correction we do not use a direct summation, but instead solve the equation 

(H - e)|*V«) = -Vi\ns), e = e ns (15) 

for each particular \ns) state. Here H is the Dirac Hamiltonian with the potential (JJJ. Similarly, the bottom left leg 
of the same diagram is evaluated using 

(H-e ns )8ij d = -V a \ns}. (16) 
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In solving this equation we take the finite size of the nucleus into account by replacing the <5-function in @ with a 
realistic nuclear density. 

Apart from the coefficients presented in Fig. which in essence show the number of diagrams of each kind, each 
particular diagram in Fig. [3] contributes with its own angular coefficient. In calculating the coefficients we assumed, 
without loss of generality, that the total angular momentum of the 4/-electrons is directed along the z-axis, \J,J Z ). 
Values of the coefficients are presented in Table IVTTT1 in the Appendix. The method for separating the radial equations 
corresponding to (|15|) and (|16[) is also described in the Appendix. As the result of the calculations we find the following 
P-odd energy correction related to the displacement Ar of the RE impurity ion: 

Ae = K a aA—(Ar -[I x J])E . (17) 
a B 

We recall that I is spin of the nucleus, J is the total angular momentum of the /-electrons, Eq — 27.2 eV is the atomic 
unit of energy, og is the Bohr radius, a is the fine structure constant, and K a is given in Eq. The dimensionlcss 
coefficient A for the Pr 3+ and Tm 3+ ions (in the corresponding lattices) calculated at D = 1.0 in Eq. (JfJJ) reads: 

A Pr = -25.99 - 11.20 + 0.32 + 0.59 + 18.64- 18.99 
+ 0.58 + 1.39 - 15.99 + 28.37 + 25.73 = 3.45, 



A Tm = 9.77+ 12.78- 3.58- 1.33 - 32.24 + 36.49 

+ 0.21 + 0.12 + 53.48 - 67.70 - 10.95 = -2.95. (18) 

The eleven terms in <|18H represent the contributions of the eleven diagrams in Fig. [21 As one can see, there is 
significant compensation between different terms in (|18fl . This compensation is partially related to the fact that each 
particular diagram in Fig. [^contains contributions forbidden by the Pauli principle. These contributions are canceled 
out only in the sum of the diagrams. To check (|18(l we have also performed a more involved calculation explicitly 
taking into account the Pauli principle in each particular diagram, the results read: 

A Pr = -0.07- 0.18 + 1.48 + 0.72 + 1.00- 2.63 

- 1.74 + 0.88 + 41.34 - 40.00 + 2.65 = 3.45, 



A Tm = -0.55- 0.56 - 4.09- 1.36 - 1.47 + 6.27 

+ 0.53 + 0.15 - 38.94 + 38.06 - 0.99 = -2.95. (19) 

Although each individual term has changed compared to l|18fl . the total sum of the diagrams remains the same. 
Comparison between l|18f) and l|19|l is a test of the many body perturbation theory used in the calculation. To 
demonstrate the sensitivity to parameters of the effective potential, we plot in Fig. 21 the coefficient A versus the 
width D of the oxygen potential, see Eq. ©. As we pointed out in Section III, the most realistic value of D is around 
0.5-1.0 . To be specific, in the final estimates we use the results (|18|) and l|19(l . which correspond to the value D = 1.0 . 




FIG. 4: Value of the coefficient A defined in Eq. I|17|l versus width of the effective oxygen potential. The dashed line corresponds 
to Pr 3+ in YGG and the solid line corresponds to Tm 3+ ions in YAG. 
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V. CALCULATION OF THE EFFECTIVE HAMILTONIAN © 



The P-odd effective Hamiltonian considered in the previous section arises due to a shift of the environment with 
respect to the rare earth ion. In other words, it is due to the first harmonic in the electron density induced by the 
perturbation operator V\ l|llf) . In the equilibrium position the first harmonic vanishes identically due to the symmetry 
of the lattice. The next harmonic in the electron density that contributes to the parity nonconserving effect is the 
third harmonic which is nonzero even in the equilibrium position of the rare earth ion. This effect gives the P-odd 
energy shift even in the absence of an external electric field. 

The effective oxygen potential Vo © represents the spherically symmetric part of the real potential for electrons 
created by the eight oxygen ions in the garnet lattice. Let us describe the potential (pseudopotential) of a single 
oxygen ion as g5(r — R), where R is the position of the ion and g is some constant. Then the total potential is 

V(r)=J29S(r-R), (20) 

R 

where summation is performed over the coordinates of the eight oxygen ions presented in Table ^ Expanding the 
Dirac delta function in the potential V(r) in a series of spherical harmonics, we find 

V(r) = g S ±^lj2J2 Y ^(R) ■ YUr). (21) 

km Ft 



Then, 



Vo(r) = 9 6(r R2 R) E Y oo(R) ■ Y oo(r) -> -A a e-^)\ (22) 

R 



and hence the third harmonic reads 



V 3 (r) = g S -^-Mj2 Y Sm(R) ■ Y 3m (r) - -A e-^) 2 ^ T 3m ■ Y 3m (r), (23) 



P 2 



R 



T 3m - y^y 3 * m (-R)- 



R 

The spherical tensor T 3m (lattice octupole) for yttrium aluminium garnet and yttrium gallium garnet has only one 
non-zero independent component, T31 = —0.1876 for YAG and T31 = —0.1010 for YGG. All other components are 
determined by the following relations: 

P33 = y^Psi, P3-1 = — Tsi, P3-3 = — P33, P30 = 0. (24) 

Components of the corresponding Cartesian irreducible tensor Tkim can be found using the following relations: 

1 g 

TxZZ T ZXZ T ZZX T X yy Ty X y Tyy X \ l ~~T 3 l . (^5) 

V J-0 

All other components of the Cartesian tensor are equal to zero. 

Similar to the "dipole" effect considered in the previous section, the octupole effect arises in the third order of 
perturbation theory. The relevant perturbation theory operators are a) interaction of the electron with the nuclear 
anapole moment V a ©, b) interaction of the electron with the lattice octupole harmonic V3 H23JI . and c) the residual 
electron-electron Coulomb interaction Vc <|12|) . The formula for the energy correction (|13fl yields 7 diagrams which 
are presented in Fig. [SJ 

Besides the coefficients presented in Fig. which show the number of diagrams of each kind, each particular 
diagram in Fig.0contributes with its own angular coefficient. In calculating the coefficients we assumed, without loss 
of generality, that the total angular momentum of 4/ electrons is directed along the z-axis, | J, J z ), and the nuclear 
spin is directed along the y-axis, / = (0,7,0). The angular coefficients for each of the 7 diagrams from Fig. [S] are 
presented in Table IVIIII in the Appendix. The method for separating the radial equations is also described in the 
Appendix. The effective Hamiltonian for the lattice octupole effect has the following form 



Ae — K a aBIi€ijkTklm {JjJiJm + JmJiJj) E . 



(26) 
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FIG. 5: Diagrams for the "octupole" effect. The cross denotes the anapole interaction V a ©, the dashed line denotes the lattice 
octupole V3 l|24|l . and the wavy line denotes the Coulomb interaction Vc The multipolarity of the Coulomb interaction 

is shown near the wavy line. Each diagram contributes with the coefficient shown before the diagram (number of diagrams of 
this kind). Summation over all intermediate states \k) and \m) and over all filled states \n) is assumed. 



Eq. (|26J) represents the only P-odd scalar combination one can construct from the two vectors and one irreducible 
third rank tensor. Note, that J here is an operator, and different components of J do not commute. This is why 
in the right hand side of Eq. I|26|) we explicitly write the Hcrmitian combination. The matrix element of (|26|l in the 
kinematics which we consider for the calculation of the angular coefficients (Table IVIHf) is 

(J, Jz\Iie ijk Tkim {JjJlJm + JmJlJj) \J, J,) = T ZZX IJ Z [5J 2 Z - 3J(J + 1) + 1] . (27) 

Our calculations show that contributions of the diagrams with the intermediate /-state (diagrams 4,5,6,7 in Fig. 03) 
are at least 30 times smaller compared to diagrams 1 and 2. The reason for this is very simple: /-electrons are 
practically decoupled from the lattice deformation. The diagram 3 is even smaller because internal 3d- and 4d- 
electrons are also decoupled from the lattice. So, only diagrams 1 and 2 contribute to the effect and they are 
nearly saturated by the intermediate unoccupied 5d-state. The dimensionless coefficient B for Pr and Tm ions in 
corresponding lattices calculated at D — 1.0 [Eq. |JBJ] reads: 

B PT (D = 1) = -2.18 + 0.76 = -1.42, 

B Tm {D = l) = 1.11 - 0.45 = 0.66. (28) 

The two terms in equations 1)28(1 represent the contributions of the first and second diagrams. The variation of the 
coefficient B with the width of the effective oxygen potential D is shown in Fig. El Again, we recall that the most 
realistic value of D is around 0.5-1.0 . To be specific, in the estimates for the effect we use D = 1.0 . 



VI. CRYSTAL FIELD, AVERAGE ELECTRON MAGNETIZATION, ORIENTATION OF NUCLEAR 

SPIN 

The energy of a free ion is degenerate with respect to the z-projection of total angular momentum. Interaction 
with the lattice (crystal field) breaks the rotational invariance and lifts the degeneracy. The effective crystal-field 
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a; 




FIG. 6: Value of the coefficient B in Eq. ® versus the width of the effective oxygen potential. Dashed line corresponds to 



Pr a+ in YGG and solid line corresponds to Tm 3+ in YAG 



Hamiltonian can be written in the following form, see, e.g., [27 



H cf = £ B* km ]T yj ^L. Ykm {r) (29) 

km i 

where Bk m are the crystal field parameters, r is the radius-vector of the atomic electron. 

Experimental values of the energy levels for Pr 3+ in YGG and Tm 3+ in YAG are known 0, 0] , and fits of the 
crystal field parameters Bkm have been performed in the experimental papers. Unfortunately, we cannot use these 
fits because they are performed without connection to a particular orientation of crystallographic axes. We need to 
know the connection and therefore we have performed independent fits. For the fits we use a modified point-charge 
model. In the simple point-charge model the crystal field is of the form 



3 3 

B ( Z ] = PkA£\ (31) 

where j enumerates ions of the lattice and pk = (r ) is the expectation value over the RE /-electron wave function. 
The values of pk are known p?T |. The point charges are qp — — 2 and qy = qca = QAl = 3. Clearly, the naive 
point-charge model is insufficient to describe the nearest 8 oxygen ions because of the relatively large size of the ions 
(extended electron density of the host oxygens) . To describe the effect of the extended electron density we introduce 
an additional field 

Akrn = A^J + A^J, (32) 



4£ - -^EJW^ y -^ ^ 



r 

3=1 3 



2n + l 



here the sum runs over the eight oxygen ions surrounding the dopant ion in the garnet structure, and are fitting 
parameters. So, we have only three fitting parameters, 0*4, and ae, because higher multipoles do not contribute 
in /-electron splitting. In the end, we get a fairly good fit of the experimental energy levels, see Table [S The values 
of the resulting crystal field parameters B^m are presented in Table IVll 

For the non-Kramers ions, such as Pr 3+ and Tm 3+ , the expectation value of the total angular momentum in the 
ground state vanishes due to the crystal field, (J) = 0. To get a nonzero (J) one needs to apply an external magnetic 
field B. Diagonalizing the Hamiltonian matrix of the dopant ion in the magnetic field 

(J' z \H cf +p B g(JB)\J z ), (34) 

(9x9 matrix for Pr 3+ and 13 x 13 matrix for Tm 3+ ) we find the ground state of the ion in the presence of the 
external magnetic field B (here pg is the Bohr magneton and g is the atomic Lande factor; g = 0.80 for Pr 3+ in 
3 H4 configuration and g = 1.17 for Tm 3+ in 3 Hq configuration.) For weak magnetic field the average total angular 
momentum can be written as 



(Jj) — TikBk- 



(35) 



11 



Pr a+ :YGG Tm 3+ :YAG 



Exp. [15] 


Calc. 


Exp. [14] 


Calc. 














23 


23 


27 


27 


23 


23 


216 


182 




400 


240 


240 


532 


413 


247 


253 


578 


538 


300 


301 


598 


621 


450 


306 


626 


877 


588 


494 


689 


895 


610 


609 






650 


673 






690 


686 






730 


825 








937 



TABLE V: Experimental and calculated crystal field energy levels in cm 1 . J- J mixing is neglected in the calculation. 



Compound 


B20 B21 


B22 


B40 


B41 


B42 


B43 


£44 


B 60 


Bei 


B 6 2 


B a3 


#64 


B 65 


B 66 


Pr:YGG 


622 Hi 


-762 


211 


-475i 


727 


1256i 


-423 


963 


-280i 


-648 


-437i 


91 


304i 


-961 


Tm:YAG 


257 92i 


-315 


-1198 


344i 


-248 


-909i 


-523 


-938 


528i 


569 


816i 


94 


-563i 


843 



TABLE VI: Crystal field parameters in cm , that fit the energy levels in Table IVI 



The tensor can be diagonalized. According to our calculations, both for Pr and Tm it is diagonal with the principal 
axes m = (1,0,0), n 2 = (1,1/^,1/^2), n 3 = (1, l/y/2, -l/y/2): 



Pr 



-0.003 
-0.154 














-0.176 



tesla 



Tm 



-0.474 
-0.023 














-0.032 



tesla 



(36) 



The average total electron angular momentum in the magnetic fields applied along the directions n 1; n 2 , and n 3 is 
plotted in Fig. [7| We see that the linear expansion (|33|) is valid for the field B < 5-10 T. 




FIG. 7: The average total electron angular momentum of the rare earth ion versus magnetic field (tesla). Directions of the 
magnetic field correspond to the principal axes of the magnetization tensor m, ri2, and n-i,. Solid lines correspond to Tm 3+ in 
YAG and dashed lines correspond to Pr in YGG. 



The effective Hamiltonian for the nuclear spin is 



A hf (J.I)-^(B-I), 



(37) 



where A^f is the hyperfine constant, fj, is the nuclear magnetic moment in nuclear magnetons and /ijy is the nuclear 
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magneton: 

141 Pr : A hf = 1093 MHz gg, fj, = 4.2754 H, / = 5/2; 
169 Tm : A hf = -393.5 MHz M = -0.2316 ||, 7=1/2. (38) 



Pr 



Tm 




9 = 0,7t 

e = 7t/4 
e = 3^/4 



9 = jc/2 




9 = 7t/2 

6 = 71/4 
6 = 37C/4 



e = o,7t 



FIG. 8: The NMR frequency versus the orientation of magnetic field with respect to the crystallographic axes, B = 10 T. We 
show the dependence on <j> for different values of 6. 

Equation 1)37(1 . together with l|35|) . gives the NMR frequency Dependence of the frequency on the orientation of 
the magnetic field B = B(s'm9 cos(j),s'm9s'm(j),cos9) with respect to the crystallographic axes is plotted in Fig. El we 
take B — 10T. Equation l|37|l also defines the quantization axis for the nuclear spin: 

n, oc (A hf (J) = (A hf fB - . (39) 

This allows us to find cross product rti x J that appears in the anapole induced energy correction (JIJ, l|17|): 



M = nj x J 



l/AftXfB) - w N B\ 



(40) 



The value of M depends on the magnitude and the orientation of the external magnetic field B with respect to the 
crystallographic axes. At B = 10 T the maximum value of M is 



Pr : M= 1.02- NT 1 , 
Tm : M = 0.79-1Q- 1 . 



(41) 



Unfortunately, the values of M are relatively small compared to the maximum possible value M = J (4 for Pr and 6 
for Tm). The suppression is due to the fact that in the nuclear magnetic Hamiltonian (|37|l the hyperfine interaction 
Ahf(J • I) is an order of magnitude larger than the direct magnetic interaction hhn(B • I)/ 1, while to maximize 
M one has to have these interactions comparable. In spite of the suppression, the observable effects related to the 
effective Hamiltonian (QJ, i|17|) are quite reasonable (see next Section). 

The situation with the effective interaction J3J|, l|2t)|) is different. Looking at equations JSJl, jUSt one can expect at 
first sight that the corresponding energy shift is nonzero only if I x (J) ^ 0. However, this is incorrect. The point 
is that due to the crystal field the tensor (JjJiJ m + J m JiJj) has nonzero components orthogonal to (J). And the 
octupole induced energy shift is in fact maximum when J || (J). The dependence of the kinematic coefficient (see 
Eq. d)J}) 



N - jli^ijk^kl'rni^J j JlJm ~\~ JjnJlJj) 



(42) 



on the orientation of magnetic field B = B(sm 9 cos (f>, sin 9 sin <f>, cos 9) at B = 10 T is plotted in Fig. [5] The maximum 
value of N is 



Pr : N = 1.81, 
Tm : N = 2.42. (43) 
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4> * 

FIG. 9: The kinematic coefficient N 1421 in the lattice octupole induced energy correction versus orientation of magnetic field 
with respect to the crystallographic axes, B = 10 T. We show the dependence on cj> for different values of 6. 

The calculations in the present section are based on the fit of experimental energy levels, Table[V] using the crystal 
field parameters. We use the set of parameters presented in Table IVll Unfortunately, the set is not unique and there 
are other sets which also reasonably fit the energy levels. In particular, for Tm 3+ in YAG there is a set of parameters 
which gives a lattice octupole induced PNC energy shift an order of magnitude larger than the present set. At this 
stage we prefer to continue with the conservative estimate. To elucidate the uncertainty related to the crystal field 
parameters detailed measurements of NMR frequencies, as well as transition amplitudes, are necessary. 



VII. ESTIMATES OF OBSERVABLE EFFECTS 



The effect ijljl. (1171) requires a displacement of the impurity ion from its equilibrium position. Such displacement can 
be achieved by application of an external electric field. The displacement has been estimated in Ref. in relation 
to the discussion of electric dipole moments. The idea behind the estimate is very simple. Since the Ga-0 link in 
YGG and the Al-0 link in YAG are much more rigid than the Y-0 links (see discussion in the electrostatic 
polarization in YGG and YAG is mainly due to displacement of the yttrium ions 

P = 3enAr. (44) 

On the other hand, the dielectric polarization caused by the external electric field E is 

P = e -=±E, (45) 

where the static dielectric constant is e « 12 for YGG and YAG. This yields the following expression for the displace- 
ment of the yttrium ions: 

e — lE 

Ar = 



4-7T 3en 



Ar/a B = 3.0 • 1(T 8 M— ]. (46) 
cm 

Measurements of infrared spectra, as well as measurements of the dependence of the dielectric constant on the 
concentration of impurities, can help to improve the estimate <|46[) . 

Using (|17|) . together with l|41|l and 1)46(1 . we obtain the following estimates for the NMR frequency shift (AI = 1) 
due to the nuclear anapole moment: 

V 

Pr : Av ~ 0.9 x 10 _9 M — 1 Hz, 

cm 

Tm : Av ~ 0.5 x 10~ 9 £[— ] Hz. (47) 



An alternative possibility for the experiment is to provide the maximum possible value of the cross product I x J 
by applying an RF pulse and then to measure the induced electric field. Using l|17|l . together with estimates of the 
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elastic constant with respect to the shift of the rare earth ion performed in [lg, we arrive at the following values of 
the anapole induced electric field: 

Pr : £~1.4x 1CT 6 V/cm, 
Tm : £~ 0.4 x 1CT 6 V/cm. (48) 

The field precesses around the direction of the magnetic field with a frequency of about 1 GHz due to the nuclear 
spin precession. In the estimates (|48|l we assume that all yttrium ions are substituted by the rare earth ions. 

Another manifestation of nuclear anapole moment is the static perpendicular macroscopic magnetization induced 
by an external electric field, 

51 <xB x E. (49) 

The exact value of the macroscopic magnetization depends on temperature and other experimental conditions, there- 
fore we cannot present a specific value. However, we can compare the effect with that expected in the electron EDM 
experiment 01 (correlation SJ cx E) using the present experimental limit on d e ^(|, d e = 1.6 x 10~ 27 e cm, as a 
reference point. The effective anapole interaction (|17|) is four order of magnitude larger than the similar effective 
EDM interaction O n the other hand, the electron EDM interaction causes electron magnetization whereas the 
anapole interaction causes only nuclear magnetization, so we lose 3 orders of magnitude on the value of the magnetic 
moment. Therefore, altogether, one should expect that the anapole magnetization is several times larger than the 
EDM magnetization. 

The effective interaction (I26[) is independent of the external electric field and is due to the asymmetric environment 
of the rare earth ion site. Since there is always another site within the unit cell which is the exact mirror reflection 
of the first one, the energy correction l|26[> does actually lead to the NMR line splitting. Using Eqs. (|26|l . (|28|l . 142fl . 
and l|43[) . we find the maximum value of this splitting corresponding to the magnetic field B = 10 T: 

Pr : Av ~ 0.5 Hz, 
Tm : Av ~ 0.25 Hz. (50) 

The splitting depends on the orientation of the magnetic field with respect to the crystallographic axes, see Fig. |SJ 



VIII. CONCLUSIONS 



In the present work we have considered effects caused by the nuclear anapole moment in thulium doped yttrium 
aluminium garnet and praseodymium doped yttrium gallium garnet. There are two effects related to the frequency 
of NMR: 1) NMR line shift in combined electric and magnetic fields. The shift is about 10~ 5 Hz at B — 10 T and 
E = 10 kV/cm. 2) NMR line splitting (magnetic field only). The spitting is about 0.5 Hz at B = 10 T. The value of 
the splitting depends on the orientation of the magnetic field with respect to the crystallographic axes. Another PNC 
effect is the induced RF electric field orthogonal to the plane of the magnetic field and nuclear spin, E oc [B x I]. 
The field is E ~ 10~ 6 V/cm at magnetic field B — 5-10 T. The last effect we have discussed is unrelated to NMR. 
This is a variation of the static macroscopic magnetization in combined electric and magnetic fields, SM oc B x E. 
The magnitude of the effect is several times larger than that expected in the electric dipole moment experiment |13| . 

It is our pleasure to acknowledge very helpful discussions with D. Budker, V.V. Yashchuk, A.O. Sushkov and A.I. 
Milstein. 



IX. APPENDIX. RADIAL EQUATIONS 



In order to calculate the energies and wavefunctions of unperturbed states of the single impurity ion in the garnet 
environment, we use the Dirac equation 

(H-e)\yj)=0. (51) 

The effective potential V(r) J7J in the Dirac Hamiltonian H is spherically symmetric, and thus the two-component 
wavefunction lib) is of the form 



if f(r)n. 

r \iag(r)fl,. 



(52) 
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Vp = Va, © 

\ip) = \ns 1/2 ) 



V p = Va,® 
\ip) = \np 1/2 ) 



= \ns 1/2 ) or 1-0) = \n Pl/2 ) 



Vp = Vs,® 

\ip) = \ns 1/2 ) or \ip) = \np 1/2 ) 



R f -K a ap n (x)f(x) 
R g \K a ap n (x)g(x) 



-\K a ap n {x)f{x) 
KaOtp„(x)g(x) 



) (r-r ) 



A a e 



fix) 



A a e (■ L rP i ) a 2 g(x) 
-Aoe-i^T f(x) 



-2i(I ■ j) 



W ■ 3) 



Ax sin 9 cos (f> + Ay sin sin (f> + Az cos 8 



jT 3m ■ Y 3m (r) 



TABLE VII: The functions Rf,R g and <£» for the different perturbation operators and different states \tp). p n is nuclear density 
normalized to unity. 



Diagram 



Pr 



Tm 



1,7,8,11 

2,9,10 

3,4,5,6 

1,2,3 

4,5 

6,7 



Dipole effect 



3^(2) - 3^LF(4) 
3^2_F(2) + J^F(4) 



3-5-7 



F(3) 



Lattice octupole effect 



T 13 ' F . (3) 

31 V 21 2-3^ -5^ -7-11 



^F(2)-^F(4) 
-3^7^(3) 



31 V 21 2^-3-5-7-H 



13-29-F(2) , 5-13-F(4) 
3^.53.72.!! + 2-3 4 -7 a -ll 2 



r 31 



13 -F(2) 



21 I 2^-3-5 3 -7 Li -ll i" 2 ;4 -3 3 -5-7 ;i -ll ;4 



13-47-F(4) 



To1 r f(2) F(4) 

J31 V 21 L2-3-5^-7-ll 2^-3 a -5-7-ll 

31 V 21 |_2-3-5^-7-ll 



2^-F(4) 
3^-5-7-11^ 



TABLE VIII: Dipole effect: Angular coefficient for each of the 11 diagrams shown in Fig. El The factor {Axl v — AyI x )J z , 
which corresponds to the kinematic structure 11711 and which is common for all the contributions, is omitted. F(l) denotes the 
Coulomb integral of multipolarity I in the radial part of the diagam. 

Lattice octupole effect: Angular coefficients for each of the 7 diagrams shown in Fig. The factor IJ Z [5J 2 - 3J(J + 1) + 1], 
which corresponds to the kinematic structure 1261 and which is common for all the contributions, is omitted. 



Here Q K and f2 K are the spherical spinors and /(r) and g(r) are radial wavefunctions. Substituting expression 1521) 
for \ip) into the Dirac equation lt5l"l) . one gets the following radial equations 

/' + K ,f/x+(-2 + a 2 (V -e))g = 0, 

g' - kq/x -{V- e)f = 0. (53) 

Here x — r/as is the radius in atomic units; k — (— iy +1 ^ 2 ~ l {j + 1/2), where j and I are the total and orbital 
angular momenta of the single-electron state correspondingly; the potential V{x) ,as well as the energy e, is expressed 
in atomic energy units. Solving the system of equations (|53|l as an eigenvalue problem numerically on a logarithmic 
coordinate grid, we find energies and wavefunctions of the unperturbed states. 
The inhomogeneous Dirac equations (|15fl and (|l(j|) are of the form 

{H-E)\5i>)=-%\^), (54) 

where V p is the single-particle perturbation operator. The correction \6if>) is of the form 



and hence the corresponding radial equations are 

F' + k'F/x + (-2 + a 2 {V - e))G = R f (n K ,\$\n K ), 

G' - k'G/x -(V- e)F = R g {n K ,\$\n K ). (56) 

The operator $ represents the angular part of the perturbation V p , and Rf and R g are the radial parts of the 
perturbation. The functions Rf, R g , and $ for all the cases we need in the present work are presented in Table IVlTl 

Having separated the radial parts, one can calculate the angular coefficients for the diagrams in Figs. Eland El The 
results of these calculations are presented in Table IVTTT1 The electronic configurations of Pr 3+ and Tm 3+ are similar: 
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two /-electrons in Pr 3+ and two /-holes in Tm 3+ . However, their orbital and spin angular momenta combine to yield 
different total angular momenta, and this makes the angular coefficients for Pr 3 + and Tm 3+ different. 
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